Human habitat modification, not apex scavenger decline, drives isotopic niche variation in a carnivore community.

Top carnivores can influence the structure of ecological communities, primarily through competition and predation; however, communities are also influenced by bottom-up forces such as anthropogenic habitat disturbance. Top carnivore declines will likely alter competitive dynamics within and amongst sympatric carnivore species. Increasing intraspecific competition is generally predicted to drive niche expansion and/or individual specialisation, while interspecific competition tends to constrain niches. Using stable isotope analysis of whiskers, we studied the effects of Tasmanian devil Sarcophilus harrisii declines upon the population- and individual-level isotopic niches of Tasmanian devils and sympatric spotted-tailed quolls Dasyurus maculatus subsp. maculatus. We investigated whether time since the onset of devil decline (a proxy for severity of decline) and landscape characteristics affected the isotopic niche breadth and overlap of devil and quoll populations. We quantified individual isotopic niche breadth for a subset of Tasmanian devils and spotted-tailed quolls and assessed whether between-site population niche variation was driven by individual-level specialisation. Tasmanian devils and spotted-tailed quolls demonstrated smaller population-level isotopic niche breadths with increasing human-modified habitat, while time since the onset of devil decline had no effect on population-level niche breadth or interspecific niche overlap. Individual isotopic niche breadths of Tasmanian devils and spotted-tailed quolls were narrower in human-modified landscapes, likely driving population isotopic niche contraction, however, the degree of individuals' specialisation relative to one another remained constant. Our results suggest that across varied landscapes, mammalian carnivore niches can be more sensitive to the bottom-up forces of anthropogenic habitat disturbance than to the top-down effects of top carnivore decline.

Complex associations between cancer progression and immune gene expression reveals early influence of transmissible cancer on Tasmanian devils.

The world's largest extant carnivorous marsupial, the Tasmanian devil, is challenged by Devil Facial Tumor Disease (DFTD), a fatal, clonally transmitted cancer. In two decades, DFTD has spread across 95% of the species distributional range. A previous study has shown that factors such as season, geographic location, and infection with DFTD can impact the expression of immune genes in Tasmanian devils. To date, no study has investigated within-individual immune gene expression changes prior to and throughout the course of DFTD infection. To explore possible changes in immune response, we investigated four locations across Tasmania that differed in DFTD exposure history, ranging between 2 and >30 years. Our study demonstrated considerable complexity in the immune responses to DFTD. The same factors (sex, age, season, location and DFTD infection) affected immune gene expression both across and within devils, although seasonal and location specific variations were diminished in DFTD affected devils. We also found that expression of both adaptive and innate immune genes starts to alter early in DFTD infection and continues to change as DFTD progresses. A novel finding was that the lower expression of immune genes MHC-II, NKG2D and CD8 may predict susceptibility to earlier DFTD infection. A case study of a single devil with regressed tumor showed opposite/contrasting immune gene expression patterns compared to the general trends observed across devils with DFTD infection. Our study highlights the complexity of DFTD's interactions with the host immune system and the need for long-term studies to fully understand how DFTD alters the evolutionary trajectory of devil immunity.

Intergenomic signatures of coevolution between Tasmanian devils and an infectious cancer.

Coevolution is common and frequently governs host-pathogen interaction outcomes. Phenotypes underlying these interactions often manifest as the combined products of the genomes of interacting species, yet traditional quantitative trait mapping approaches ignore these intergenomic interactions. Devil facial tumor disease (DFTD), an infectious cancer afflicting Tasmanian devils (Sarcophilus harrisii), has decimated devil populations due to universal host susceptibility and a fatality rate approaching 100%. Here, we used a recently developed joint genome-wide association study (i.e., co-GWAS) approach, 15 y of mark-recapture data, and 960 genomes to identify intergenomic signatures of coevolution between devils and DFTD. Using a traditional GWA approach, we found that both devil and DFTD genomes explained a substantial proportion of variance in how quickly susceptible devils became infected, although genomic architectures differed across devils and DFTD; the devil genome had fewer loci of large effect whereas the DFTD genome had a more polygenic architecture. Using a co-GWA approach, devil-DFTD intergenomic interactions explained ~3× more variation in how quickly susceptible devils became infected than either genome alone, and the top genotype-by-genotype interactions were significantly enriched for cancer genes and signatures of selection. A devil regulatory mutation was associated with differential expression of a candidate cancer gene and showed putative allele matching effects with two DFTD coding sequence variants. Our results highlight the need to account for intergenomic interactions when investigating host-pathogen (co)evolution and emphasize the importance of such interactions when considering devil management strategies.

Disease-driven top predator decline affects mesopredator population genomic structure.

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection.

BACKGROUND: Transmissible cancers lie at the intersection of oncology and infectious disease, two traditionally divergent fields for which gene expression studies are particularly useful for identifying the molecular basis of phenotypic variation. In oncology, transcriptomics studies, which characterize the expression of thousands of genes, have identified processes leading to heterogeneity in cancer phenotypes and individual prognoses. More generally, transcriptomics studies of infectious diseases characterize interactions between host, pathogen, and environment to better predict population-level outcomes. Tasmanian devils have been impacted dramatically by a transmissible cancer (devil facial tumor disease; DFTD) that has led to widespread population declines. Despite initial predictions of extinction, populations have persisted at low levels, due in part to heterogeneity in host responses, particularly between sexes. However, the processes underlying this variation remain unknown. RESULTS: We sequenced transcriptomes from healthy and DFTD-infected devils, as well as DFTD tumors, to characterize host responses to DFTD infection, identify differing host-tumor molecular interactions between sexes, and investigate the extent to which tumor gene expression varies among host populations. We found minimal variation in gene expression of devil lip tissues, either with respect to DFTD infection status or sex. However, 4088 genes were differentially expressed in tumors among our sampling localities. Pathways that were up- or downregulated in DFTD tumors relative to normal tissues exhibited the same patterns of expression with greater intensity in tumors from localities that experienced DFTD for longer. No mRNA sequence variants were associated with expression variation. CONCLUSIONS: Expression variation among localities may reflect morphological differences in tumors that alter ratios of normal-to-tumor cells within biopsies. Phenotypic variation in tumors may arise from environmental variation or differences in host immune response that were undetectable in lip biopsies, potentially reflecting variation in host-tumor coevolutionary relationships among sites that differ in the time since DFTD arrival.

Isotopic niche variation in Tasmanian devils Sarcophilus harrisii with progression of devil facial tumor disease.

Devil facial tumor disease (DFTD) is a transmissible cancer affecting Tasmanian devils Sarcophilus harrisii. The disease has caused severe population declines and is associated with demographic and behavioral changes, including earlier breeding, younger age structures, and reduced dispersal and social interactions. Devils are generally solitary, but social encounters are commonplace when feeding upon large carcasses. DFTD tumors can disfigure the jaw and mouth and so diseased individuals might alter their diets to enable ingestion of alternative foods, to avoid conspecific interactions, or to reduce competition. Using stable isotope analysis (δ13C and δ15N) of whiskers, we tested whether DFTD progression, measured as tumor volume, affected the isotope ratios and isotopic niches of 94 infected Tasmanian devils from six sites in Tasmania, comprising four eucalypt plantations, an area of smallholdings and a national park. Then, using tissue from 10 devils sampled before and after detection of tumors and 8 devils where no tumors were detected, we examined whether mean and standard deviation of δ13C and δ15N of the same individuals changed between healthy and diseased states. δ13C and δ15N values were generally not related to tumor volume in infected devils, though at one site, Freycinet National Park, δ15N values increased significantly as tumor volume increased. Infection with DFTD was not associated with significant changes in the mean or standard deviation of δ13C and δ15N values in individual devils sampled before and after detection of tumors. Our analysis suggests that devils tend to maintain their isotopic niche in the face of DFTD infection and progression, except where ecological conditions facilitate a shift in diets and feeding behaviors, demonstrating that ecological context, alongside disease severity, can modulate the behavioral responses of Tasmanian devils to DFTD.

Quantifying 25 years of disease-caused declines in Tasmanian devil populations: host density drives spatial pathogen spread.

Infectious diseases are strong drivers of wildlife population dynamics, however, empirical analyses from the early stages of pathogen emergence are rare. Tasmanian devil facial tumour disease (DFTD), discovered in 1996, provides the opportunity to study an epizootic from its inception. We use a pattern-oriented diffusion simulation to model the spatial spread of DFTD across the species' range and quantify population effects by jointly modelling multiple streams of data spanning 35 years. We estimate the wild devil population peaked at 53 000 in 1996, less than half of previous estimates. DFTD spread rapidly through high-density areas, with spread velocity slowing in areas of low host densities. By 2020, DFTD occupied >90% of the species' range, causing 82% declines in local densities and reducing the total population to 16 900. Encouragingly, our model forecasts the population decline should level-off within the next decade, supporting conservation management focused on facilitating evolution of resistance and tolerance.

A transmissible cancer shifts from emergence to endemism in Tasmanian devils.

Emerging infectious diseases pose one of the greatest threats to human health and biodiversity. Phylodynamics is often used to infer epidemiological parameters essential for guiding intervention strategies for human viruses such as severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2). Here, we applied phylodynamics to elucidate the epidemiological dynamics of Tasmanian devil facial tumor disease (DFTD), a fatal, transmissible cancer with a genome thousands of times larger than that of any virus. Despite prior predictions of devil extinction, transmission rates have declined precipitously from ~3.5 secondary infections per infected individual to ~1 at present. Thus, DFTD appears to be transitioning from emergence to endemism, lending hope for the continued survival of the endangered Tasmanian devil. More generally, our study demonstrates a new phylodynamic analytical framework that can be applied to virtually any pathogen.

Infectious disease and sickness behaviour: tumour progression affects interaction patterns and social network structure in wild Tasmanian devils.

Infectious diseases, including transmissible cancers, can have a broad range of impacts on host behaviour, particularly in the latter stages of disease progression. However, the difficulty of early diagnoses makes the study of behavioural influences of disease in wild animals a challenging task. Tasmanian devils (Sarcophilus harrisii) are affected by a transmissible cancer, devil facial tumour disease (DFTD), in which tumours are externally visible as they progress. Using telemetry and mark-recapture datasets, we quantify the impacts of cancer progression on the behaviour of wild devils by assessing how interaction patterns within the social network of a population change with increasing tumour load. The progression of DFTD negatively influences devils' likelihood of interaction within their network. Infected devils were more active within their network late in the mating season, a pattern with repercussions for DFTD transmission. Our study provides a rare opportunity to quantify and understand the behavioural feedbacks of disease in wildlife and how they may affect transmission and population dynamics in general.

Ecological and Evolutionary Consequences of Anticancer Adaptations.

Cellular cheating leading to cancers exists in all branches of multicellular life, favoring the evolution of adaptations to avoid or suppress malignant progression, and/or to alleviate its fitness consequences. Ecologists have until recently largely neglected the importance of cancer cells for animal ecology, presumably because they did not consider either the potential ecological or evolutionary consequences of anticancer adaptations. Here, we review the diverse ways in which the evolution of anticancer adaptations has significantly constrained several aspects of the evolutionary ecology of multicellular organisms at the cell, individual, population, species, and ecosystem levels and suggest some avenues for future research.

Blood Parasites in Endangered Wildlife-Trypanosomes Discovered During a Survey of Haemoprotozoa from the Tasmanian Devil.

The impact of emerging infectious diseases is increasingly recognised as a major threat to wildlife. Wild populations of the endangered Tasmanian devil, Sarcophilus harrisii, are experiencing devastating losses from a novel transmissible cancer, devil facial tumour disease (DFTD); however, despite the rapid decline of this species, there is currently no information on the presence of haemoprotozoan parasites. In the present study, 95 Tasmanian devil blood samples were collected from four populations in Tasmania, Australia, which underwent molecular screening to detect four major groups of haemoprotozoa: (i) trypanosomes, (ii) piroplasms, (iii) Hepatozoon, and (iv) haemosporidia. Sequence results revealed Trypanosoma infections in 32/95 individuals. Trypanosoma copemani was identified in 10 Tasmanian devils from three sites and a second Trypanosoma sp. was identified in 22 individuals that were grouped within the poorly described T. cyclops clade. A single blood sample was positive for Babesia sp., which most closely matched Babesia lohae. No other blood protozoan parasite DNA was detected. This study provides the first insight into haemoprotozoa from the Tasmanian devil and the first identification of Trypanosoma and Babesia in this carnivorous marsupial.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour.

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape.

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Phylogeographic structure across one of the largest intact tropical savannahs: Molecular and morphological analysis of Australia's iconic frilled lizard Chlamydosaurus kingii.

The spectacular threat display of the savannah specialist Australo-Papuan frilled lizards has made them one of the world's most iconic reptiles. They are increasingly used as a model system for research in evolutionary biology and ecology but little is known of their population structure. Their distribution across northern Australia and southern New Guinea also provides an opportunity to examine biogeographic patterns as they relate to the large-scale movement of savannah habitat during the Plio/Pleistocene and the associated increase in aridity. We generated sequence data for one mitochondrial and four nuclear DNA loci (5052 base pairs) for 83 frilled lizards sampled throughout their range. We also quantified body proportion variation for 279 individuals. Phylogenetic analyses based on maximum likelihood and Bayesian species-tree methods revealed three shallow clades that replace each other across the monsoon tropics. We found the expected pattern of male biased sexual size dimorphism in both maximum body size and head size but there was no sexual dimorphism in overall body shape or in frill size, relative to head size, supporting the hypothesis that the frill is used primarily as a threat display rather than a sexual display. The genetic clades are broadly consistent with known clinal variation in frill color that gradually shifts from west to east (red, orange, yellow/white) but otherwise show little morphological differentiation in body proportion measures. The biogeographic breaks between clades occur at the Carpentaria Gap and the lowlands surrounding the Ord River, and our ecological niche modeling predicts lower habitat suitability for C. kingii in these regions. While this biogeographic pattern is consistent with numerous other taxonomic groups in northern Australia, the overall low genetic diversity in frilled lizards across the entire monsoon tropics and southern New Guinea contrasts starkly to patterns seen in other terrestrial vertebrates. Extremely low intra-clade genetic diversity over vast geographic areas is indicative of recent gene flow that would likely have been facilitated by widespread savannah during interglacials, or alternatively may reflect population bottlenecks induced by extreme aridity during Pleistocene glacials. The shallow divergence between Australian and New Guinean samples is consistent with recent connectivity between Australia and New Guinea that would have been possible via a savannah corridor across the Torres Strait. Based on our molecular and morphological data, we do not support taxonomic recognition of any of the frilled lizard clades and instead consider C. kingii a single species with shallow phylogeographic structure and clinal variation in frill color.